Avalanches on a Conical Bead Pile: Scaling with Tuning Parameters 
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Uniform spherical beads were used to explore the scaling behavior of a granular system near its critical an- 
gle of repose on a conical bead pile. We found two tuning parameters that could take the system to a critical 
point. The existence of those tuning parameters violates the fundamental assumption of self-organized critical- 
ity, which proposed that complex dynamical systems self-organize to a critical point without need for tuning. 
Our avalanche size distributions were well described by a simple power-law with the power T = 1.5 when drop- 
ping beads slowly onto the apex of a bead pile from a small height. However, we could also move the system 
from the critical point using either of two tuning parameters: the height from which the beads fell onto the top 
of the pile or the region over which the beads struck the pile. As the drop height increased, the system did not 
reach the critical point yet the resulting distributions were independent of the bead mass, coefficient of friction, 
or coefficient of restitution. All our apex-dropping distributions for any type of bead (glass, stainless steel, 
zirconium) showed universality by scaling onto a common curve with T = 1.5 and a = 1.0, where 1/a is the 
power of the tuning parameter. From independent calculations using the moments of the distribution, we find 
values for T = 1.6±0. 1 and cr = 0.91 ±0.15. When beads were dropped across the surface of the pile instead of 
solely on the apex, then the system also moved from the critical point and again the avalanche size distributions 
fell on a common curve when scaled similarly using the same values of X and <7. We also observed that an hep 
structure on the base of the pile caused an emergent structure in the pile that had six faces with some fee or hep 
structure. 

Keywords: bead pile; avalanches; scaling behavior; granular material; universality; critical behavior; critical exponents; self- 
organized criticality 



INTRODUCTION 



Some have proposed that granular (and other) systems self- 
organize to true critical points 01 0] like dry sand forming a 
conical pile at its critical angle of repose; others have sug- 
gested that some tuning parameters must exist in order to 
bring a granular system to the critical point |3j. In this pa- 
per, we present two such tuning parameters and the remark- 
able scaling behavior associated with them. Critical points 
are characterized by a simple power-law behavior for relevant 
properties when brought sufficiently near the critical point by 
varying some tuning parameter like temperature or concentra- 
tion. In many systems (liquids, ferromagnets, polymer solu- 
tions, etc.), a wide range of scaling behavior occurs as char- 
acterized by Widom ]4|] and later formalized by renormaliza- 
tion group theory pO. The particulars of the system determine 
where the critical point is located, but when sufficiently near 
the critical point then the system's properties (like correlation 
length or susceptibility) scale with temperature as a power-law 
with powers that are universal values for all types of systems 

A granular system behaves in some ways like a liquid with 
an ability to flow and in some ways like a solid with a sta- 
ble fixed structure if undisturbed. Bak, Tang, and Wiesen- 
feld |H Q] proposed that granular piles would self-organize 
unaided by any tuning parameter to a critical point and would 
avalanche with size and temporal distributions that would be 
described by simple power-laws. As grains are added to a pile, 
a much larger percentage of avalanches would be small rather 
than large and the avalanche size distribution would follow a 



power-law relation 



P(s)=Pof 



(1) 
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where P(s) is the fraction of avalanches of size s, P () is a con- 
stant, and T is the power. The power is predicted by mean- 
field theory (MFT) to be 3/2 when the dimension of the sys- 
tem is greater than the critical dimension of four ifioll but to be 
1.33 in three dimensions fiUl . This theory of self-organized 
criticality (SOC) has been applied to many s yste ms includ- 
ing earthquakes |l 
market trends 1L5J 

The conditions under which a system exhibits SOC and 
simple power law distributions for the size and duration of its 
response have been the subject of numerous theoretical and 
numerical simulations with relatively few experimental inves- 
tigations! Jjl-Iil]; a history we have described previously in 
detail |fl9IK 

Here we have undertaken a systematic investigation of a 
simple experimental system that models a sandpile — a pile of 
spherical beads — in order to explore when the system can be 
described by a pure power-law. We have experimentally found 
two parameters that tune the system away from the critical 
point. The existence of these tuning parameters that bring a 
system to its critical point violates a fundamental assumption 
of SOC. 

In the experiment, we formed a conical bead pile on a circu- 
lar base, then dropped one bead at a time onto the pile (either 
at the apex or across the pile) and measured the number of 
beads falling off the pile, which we took to be proportional to 
the size of the avalanche. We varied the height at which beads 
were dropped, the location where they landed on the pile, the 
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type of bead, and the underlying base geometry of the pile. 
We use the analytic tools developed for equilibrium systems 
near a critical point and apply them to this non-equilibrium, 
granular system. In addition to looking at values of the criti- 
cal exponent T as the slope of a line on a log-log plot, we also 
look for universal scaling functions that can collapse our data 
onto a common curve and provide a more compelling case for 
an exponent's value. We show that these scaling functions col- 
lapse our data on several different types of beads remarkably 
well over several decades of avalanche size. 



THEORY 

As a granular system moves from the critical point and a 
pure power-law, the avalanche size distribution is typically 
characterized by a function determined from MFT in which 
the system departs from a power law at a characteristic ava- 
lanche size s a fiSHS* so that 



P(s)=P oS - T e 



(2) 



When s (1 is large, the function collapses to a simple power 
law. We have previously experimentally verified this behav- 
ior in uniform glass beads ) 1911 and found a proportional rela- 
tionship between 1 js and the height /; from which beads are 
dropped onto the apex of a near-critical bead pile. To empha- 
size the role of h as a tuning parameter, the argument of the 
exponential —s/s may be rewritten as — Ash l l a , where A is 
a proportionality constant and a is a universal exponent 121 1. 
More generally, the exponential function of the product sh}l a 
may be replaced by an unspecified scaling function 3J {shAl a \ 
for any tuning parameter h, so that the distribution becomes 



P{s)=P s~ x< &[sh 



A/a 



(3) 



A random-field Ising Model has shown the same universal 
critical exponents and scaling functions regardless of whether 
a system is in equilibrium or non-equilibrium indicating that 
these both bel ong to the same universality class in 3D 12211 . 
Sethna et al. 12311 have used a renormalization-group frame- 
work to analyze the avalanche behavior of crackling noise (in- 
cluding earthquakes, crumpled paper, and Barkhausen noise 
in magnets) and to search for the universality class govern- 
ing crackling. They emphasize that scaling functions can be 
more sensitive than critical exponents to discriminate between 
different universality classes because scaling functions con- 
tain much more information than do the small set of numbers 
which are the critical exponents. This framework provides 
specific predictions for the moments (s) and (s 2 ) of the ava- 
lanche size distribution. For a distribution given by a power- 
law modified by a scaling function (sh l / a ^ as in Eq. (01, 
then the moments of the distribution (s) and (s 2 ), defined be- 
low, may be rewritten using the scaling function and using a 



change of variables in the integration to give 
(j) = / sP(s)ds = f sP B s~ T 3i (sh 1 / ) ds oc h^, and 

Jails Alls V / 

(4a) 

- 2x - I" s 2 P(s)ds= f s 2 P s- x ^.{sh l l a ^dsoch^ 

(4b) 
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for any scaling function 3i. Thus, if the tuning parameter and 
scaling are valid, then values for the universal exponents T and 
a can be determined from the dependence of the moments of 
the distribution (s) and (s 2 ) on the tuning parameter h. 

The scaling function also provides a means of collapsing 
the data obtained at different values of the tuning parameter 
onto a single curve. Multiplying Eq. (0 by hr z l a yields 

P(s)hT z l a =P (sh l l a \ X ^ (sh 1/a ) , (5) 



so that a plot of P(s)h~ x l a versus sh l l a should collapse all the 
distribution data from different heights and bead types onto 
the same curve. Alternatively, Eq. (0 may be multiplied by 
s x instead to obtain 



P(s)s r =P tJ 3i (sh l ' a 



(6) 



which will provide a collapse on a plot of P(s)s x versus sh 1 / a . 
The two methods of collapse are complementary since we are 
interested in the values of two independent universal expo- 
nents, t and a. 



EXPERIMENTAL DETAILS 

Our experiment involves adding individual beads to a pre- 
formed conical pile of the same beads on a circular base of 
fixed diameter (17.8 cm or about 60 bead diameters). Each 
base had (with one exception described below) a monolayer of 
randomly placed beads glued to it that prevented the beads in 
the pile from rolling off the base. The pile is initially formed 
by slowly pouring beads onto the base after which beads are 
added one at a time with the equilibrium mass of the pile cap- 
tured by a Lab VIEW data acquisition program. The mass of 
the pile is measured with a precision of 0.01 g, which is suffi- 
cient to resolve bead additions or the loss of one or more beads 
(as in an avalanche). The data acquisition program waits for 
the pile mass to stabilize before recording the mass and then 
dropping the next bead. Thus the system has the very slow 
driving rate needed for SOC. The apparatus rests on a vibra- 
tion isolated, optical table and is surrounded by a plexiglass 
box to minimize air drafts and acoustic noise. In the first set 
of experiments, beads were dropped onto the apex of the pile 
from a height h that can be adjusted from 0.6 cm to 1 1 cm and 
was measured by a cathetometer with a reproducibility of 1 
mm. 

Each data run consisted of dropping approximately 20 000 
beads over a 54-hour period with multiple, reproducible runs 



3 



TABLE I: Spherical beads of different materials as used in our experiments, where fl^ and fl s are the coefficients of kinetic or static friction 
and e is the coefficient of restitution 



Material 


Diameter (mm) 


Mass (g) 


Density 
(g/cm 3 ) 


ft 


Ms 


f 


Glass (soda-lime) 


3.0 ±0.1 


0.035± 0.001 


2.5 


0.126 ± 0.005 


0.144 ±0.008 


0.90 ± 0.03 


Zirconium 


3.000 ± 0.002 


0.085 ± 0.0005 


4.49 


0.124 ± 0.005 


0.152 ±0.007 


0.76 ± 0.04 


Stainless Steel 


3.171 ±0.002 


0.133 ±0.0005 


7.28 


0.159 ± 0.008 


0.20 ± 0.01 


0.78 ± 0.04 


Steel Shot 


3.30 ± 0.04 


0.150 ±0.005 


7.28 


0.205 ± 0.006 


0.22 ± 0.01 


0.67 ± 0.05 



to allow calculation of statistical error bars. Fig. 02a) shows 
an example run of beads dropping onto the apex of the pile 
from a height of 1.5 cm. Decreases in mass correspond to 
avalanches. A narrow slice of data corresponding to 80 bead 
drops is shown enlarged in the inset; a pile-building period in- 
cluding a few one-bead avalanches and a three-bead avalanche 
can be seen at the left of the inset, followed by a large drop 
corresponding to a 30-bead avalanche, after which the pile 
again builds. An analysis routine is used to determine ava- 
lanche size and to count the number of avalanches of a given 
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FIG. 1: The mass of a conical pile of zirconium beads after each 
bead drop, (a) The beads are dropped from a height of 1.5 cm onto 
the apex, (b) The beads are dropped from a height of 2 cm (as mea- 
sured from the apex of the pile) and over a width of 7.6 cm across the 
surface. Both plots assign the vertical scale's zero arbitrarily. The in- 
set graphs expand a narrow sliver (80 bead drops) of the main graph; 
insets in both (a) and (b) have a vertical full scale of 4 g 



size and calculate the probability. (Note that one-bead ava- 
lanches are not included in our calculations of the distribution 
P(s) because they may be caused by the dropped bead falling 
immediately off the pile, by the dropped bead knocking a sin- 
gle other bead off the pile, or by a jam in the bead dropper so 
that no bead was dropped.) The typical size of our 3D bead 
pile was around 15 000 beads. Details of the equipment and 
procedure were described previously lfl9tl where we reported 
the results of a similar experiment using 3.0 mm glass spheres. 
We have repeated the experiment using beads of different den- 
sities (steel 1 and zirconium 2 ) but the same nominal diameter 
of 3 mm. Table U gives the relevant details for each of the 
beads we have used. We include our measured values for the 
coefficients of kinetic and static friction, fi^ and jj, s , as well as 
the coefficient of restitution e, which measures the inelasticity 
of the material when involved in a collision [24]. 



As mentioned above, we take the change in mass of the 
pile as proportional to avalanche size. While experiments on 
rice 12811 observed internal avalanches which did not change 
the pile mass, spherical beads have a different dynamical be- 
havior. When studying spherical glass beads on an inclined 
surface, Daerr and Douady 1E9I1 found that avalanches on the 
surface form a front that moves at an approximately constant 
velocity and incorporates the first couple of layers on the sur- 
face. Borzsonyi et al. J30I1 found a different type of flow for 
irregular material than for spherical glass beads; the latter had 
fronts with a continuous structure traveling with a front speed 
that was faster than the moving beads. Both studies involved 
beads on an incline past the critical angle of repose and, there- 
fore, indicate that avalanches in spherical beads, once formed, 
continue moving down the surface of the pile and off the edge. 

We have here investigated three experimental variations 
where identical beads were added to a pile to perturb its equi- 
librium. The first variation probed the effect of the bead's 
density on the avalanche size distribution as the drop height h 
onto the apex of the pile was varied. In the second variation, 
we explored the effect of adding beads across the surface of 
the pile. Finally, we note the emergence of structure within 
the pile when the base is a 2D hexagonal close-pack instead 
of random. 



1 New England Miniature Ball Company 

2 K-Style Advanced Ceramics, Qingdao, People's Republic of China 
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Apex dropping 

Individual beads were dropped onto the apex of the conical 
pile from an adjustable height h above the apex. We used 
both stainless steel and then zirconium beads and varied the 
drop height from 0.6 cm to 1 1 cm and from 1.5 cm to 10 cm 
respectively. Prior to use, the beads were cleaned using Micro- 
Soap and distilled water and dried thoroughly. 



large avalanches when the pile came in contact with the three 
constraining sides. Their pile then avalanched along one di- 
rection (toward the open end) and became essentially a ID 
flow while ours remained a 2D flow as the avalanches typi- 
cally occurred along the six edges on the surface of the pile. 
We did not observe quasi-periodic large avalanches in our 
experiment and the avalanche size distributions when cliffs 
formed are consistent with, and included among, the stainless 
steel data taken using a randomly packed base shown in Fig.|2] 



Pseudo-random dropping 



RESULTS 



Another experimental variation explored dropping zirco- 
nium beads, not just onto the apex, but at varied locations 
across the pile. We modified our bead dropper by adding a 
randomizer that consisted of four horizontal diamond-pattern 
metal screens vertically separated by 3.2 cm. Zirconium beads 
dropped into this randomizer bounced chaotically before pass- 
ing through the diamond-shaped openings to the next screen 
and finally emerging from the bottom at a semi-random loca- 
tion and with no vertical velocity. We measured the distribu- 
tion of beads coming out of the randomizer and found a broad 
Gaussian distribution centered above the bead pile's apex. The 
diameter of the drop area was varied from 3.2 cm up to a max- 
imum of 15.2 cm, which is slightly less than the overall pile 
diameter of 17.8 cm. The bottom of the randomizer was kept 
at a constant height of 2 cm above the apex of the pile. Beads 
falling onto the edge of the pile fell a slightly longer distance, 
up to a maximum of 5.5 cm at the edge of the largest drop 
area. 



Emergence of structure 

For one set of measurements with very uniform stainless 
steel beads, the base layer of beads was carefully glued in 
a hexagonal close-pack structure. Stainless steel beads were 
then poured onto the pile forming the typical random packing 
that resulted in an angle of repose around 22°. However, after 
many thousands of beads had been dropped onto the apex of 
the pile, we noted that a hexagonal structure emerged within 
the pile. A close-pack triangular structure had previously been 
reported in 2D piles [26] but not in 3D piles. In our piles, the 
emergent structure appeared at the bottom edge of the pile and 
worked up into a cliff at 52°or 66°that corresponded on alter- 
nating sides to face-centered cubic (fee) and hexagonal close 
packing (hep) structures. Each cliff became 3-7 bead layers 
high and all six sides associated with the hexagonal bead pat- 
tern on the base could be observed with some random packing 
above this structure. 

Interestingly, the avalanche size distribution did not change 
as a result and we observed the avalanches to occur most often 
along the six boundaries between the cliffs, which themselves 
were quite robust once they formed. This is in contrast to an 
experiment by Lorincz and Wijngaarden [27] who built a pile 
in a square box with one open side and found quasi-periodic 



Apex dropping 

For small drop heights with either bead material, the dis- 
tribution of avalanche sizes are described by the pure power 
law of Eq. (Q~|) with an exponent T of 1.5, which was the value 
found using glass beads 11911 . For these small drop heights, the 
log-log plots of the distributions were approximately linear for 
two decades in avalanche size, more than similar experiments 
but less than the three or more decades needed to determine 
an exponent precisely. As h increased, the probability P{s) 
decreased for larger avalanches. This is shown in Fig. |2a) 
for stainless steel and[2|b) for zirconium. A weighted fit to 
Eq. © can nicely describe the avalanche size distribution for 
beads dropped from larger heights as shown in Fig. |2jb). As 
with glass beads, the cutoff parameter decreases as h increases 
since s () roughly corresponds to the value of s where the prob- 
ability breaks from a pure power law. Remarkably, within our 
resolution, 1 / s c , is proportional only to h and does not depend 
on the mass of the bead (which differed by almost a factor of 
four), the coefficient of friction, or the coefficient of restitu- 
tion. 



Pseudo-random dropping 

The shape of the pile and the kind of avalanching was re- 
markably different as a result of using the randomizer. Our 
pile was initially conical from simply pouring beads onto the 
base to form the initial pile, but as beads were added by the 
randomizer the pile became flatter and had many small ava- 
lanches and far fewer larger ones. The raw mass versus time 
plot for one of these runs is shown in Fig. 02b) where the small 
variation in pile mass and the lack of clear building periods re- 
flect the many small avalanches that occurred. The change in 
pile-building behavior can be seen in more detail in the inset 
of Fig. Q2b) (set at the same scale as in Fig.QJa)); in addition 
to the few-bead avalanches visible, a substantial number of 
one -bead avalanches prevent the pile from building smoothly. 

The lack of large avalanches is apparent in the plot of the 
avalanche size distributions shown in Fig. [3] The dramatic 
roll-off of the probability was initially surprising so we var- 
ied the diameter w of the randomizer by factors of two from 
15.2 to 3.8 cm to look for systematic effects. As expected, 
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FIG. 2: The probability P(s) of avalanches of size s for (a) stainless 
steel beads and (b) zirconium beads. The probability changes from 
a pure power law for small drop heights h to being described by Eq. 
$2% for larger drop heights. Lines indicate the result of weighted fits 
to Eq. (O with z= 1.5; fit lines are not shown in (a) for clarity. Error 
bars are one standard deviation estimates 
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FIG. 3: Zirconium beads dropped from a height of 2 cm (as measured 
from the apex of the pile) but distributed over a width w of the pile. 
Apex dropping data and dashed line corresponding to a power law 
with T = 1.5 are shown for comparison. Other lines are guides for 
the eye 



as the diameter of the region over which the beads were drop- 
ped approached apex dropping, the avalanche size distribution 
approached the distribution seen for apex dropping from the 
same height. The avalanche size distribution for the smallest 
drop area, w = 3.8 cm, was well described by Eq. (O, but 
the distributions for larger w values were not well fit by this 
modified power law. 

The change in pile-building behavior we observed is some- 
what similar to that seen in a previous granular simula- 
tion which briefly investigated the effects of adding particles 
across the surface of a pile (25J. Overall the simulation did not 
produce a realistic distribution of avalanche sizes, but the au- 
thors did note that random sprinkling caused pile growth to be 
interrupted by small avalanches while central (apex) fueling 
made the pile grow uninterrupted. 



DISCUSSION 

Our bead piles reach a critical state when the added beads 
arrive with little velocity and settle onto the pile, building it 
to a locally critical state which can then avalanche with the 
addition of a single bead that subsequently carries beads off 
the surface of the pile. For beads dropped from small heights, 
we saw the mass of the pile generally increasing along with 
many small avalanches as well as the occasional large ava- 
lanche as can be seen in Fig. [TJ a )- The resulting avalanche 
size distribution was well described by a power-law. However, 
when beads were added to the apex of the pile with larger ve- 
locities, some beads that were locally sub-critical on the pile 
were dislodged resulting in more frequent small and interme- 
diate sized avalanches. These prevented the pile from build- 
ing enough for large avalanches to occur. The system was 
kept sub-critical and the avalanche size distribution rolls off 
for large avalanches as seen in Fig. [2] 

Scaling Analysis A better way of representing the behav- 
ior of the granular pile when dropping beads onto the apex is 
by scaling the distribution data as described in Sect. „ Specif- 
ically, if h is a tuning parameter, then Eq. © and Eq. © can 
be used to provide a collapse of all the distribution data at all 
drop heights and all bead types onto a common curve. This 
is shown in Fig. @] using T = 1.5 and a = 1, which best col- 
lapsed the data onto the two common curves shown in (a) and 
(b). That the avalanche size distribution for beads of very dif- 
ferent densities, coefficients of restitution, and coefficients of 
friction would collapse onto a common curve using such sim- 
ple scaling indicates that the instability structure of the conical 
pile in 3D is effectively independent of the beads used. The 
scaling parameter h allows the system to be brought close to 
the critical point when h is small but smoothly taken away 
from critical as h increases. 

A similar effect is seen when zirconium beads were drop- 
ped over a region of width w onto the bead pile from a fixed 
height. The same kind of scaling described above (when h 
was used as a parameter in Eq. (0) also works when the tun- 
ing parameter becomes w. The avalanche size distributions 
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FIG. 4: The avalanche size distribution data for all types of beads 
falling from all drop heights h onto the apex of a pile collapses onto 
a common curve when scaled using Eq. l[3} as described in the text. A 
plot of P(s)h~ x l a versus sh^l a is shown in (a); P(s)s x versus sh^l a 
is plotted in (b). The values of % = 1.5 and a = 1 were found to give 
the best collapse in both plots 
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FIG. 5: The avalanche size distributions for beads falling across the 
surface of the pile shown in Fig.[3]can be scaled according to Eq. (|5j 
(using vv in place of h) so all the data fall onto a common curve. The 
exponent values are the same as in Fig. [4] The apex data were scaled 
using w = 0.5 cm 



h and Fig. |6tb) against w. In both cases, the data show very 
good consistency with a power-law using a power of - x ji for 
(s) and ~ 3 /2 for {s 2 } as shown by the solid lines. The only 
deviation appears in the (s 2 ) values for the three smallest drop 
heights with one type (stainless steel) of bead; if those three 
points are excluded from a weighted fit to the other (s 2 ) points 
in (a) then we obtain a power of 1.55 ±0.13. A weighted fit 
to the (s) points from (a) yields a power of —0.45 ±0.05. In- 
verting these two results for z and a, we find T = 1.6 ±0.1 
and <7 = 0.91 ±0.15. This calculation is independent from 
the collapse technique of Eq. (O and Eq. © but yields values 
of exponents consistent with our finding of T = 1.5 and (7=1 
from Fig. [4] 



CONCLUSION 



shown in Fig. [5] scale onto a common curve as seen in Fig. [5] 
using this scaling approach when T = 1.5 and (7=1. This in- 
dicates that varying the region over which beads are dropped 
acts as another tuning parameter for the system in a way that 
is analogous, but not equivalent, to varying the height when 
dropping onto the apex. 

We can further test the robustness of the scaling when h or w 
act as a tuning parameter as well as testing the value of the ex- 
ponents T and <7 by calculating the two moments (s) and (s 2 ) 
of the avalanche size distribution as detailed in the literature 



|22|,|23|]. If scaling holds as discussed in Eq. (0]i then the mo- 
ments of the distribution (s) and (s 2 ) should vary as 
and h^ T ~^l a respectively, which, for T = 1.5 and <7 = 1.0, 
correspond to h~ l l 2 and hr^l 2 (or to w -1 / 2 and w~ 3 / 2 ). Using 
the discrete values from our measured probabilities, we can 
calculate the moments from each of the data runs presented 
in Figs. |2] and [3j these values are plotted in Fig. |6|a) against 



With our randomly packed, 3D conical bead piles of spher- 
ical beads on a base of approximately 60 bead diameters, 
we find the avalanche size distributions to indicate critical 
point behavior in this non-equilibrium system. At small drop 
heights onto the apex of the pile, the distributions for glass, 
zirconium, or stainless steel beads are described by a simple 
power-law with exponent T = 1.5. We find the drop height 
h and the diameter of the drop area w to be two independent 
tuning parameters that take the system away from the criti- 
cal point. Despite differences in the inertia, surface proper- 
ties, and elasticity of the beads, all the distributions collapse 
onto a common curve when scaled appropriately by the tun- 
ing parameter and using T = 1.5 and (7=1. The moments 
of those distributions show a power-law behavior in the tun- 
ing parameters that further show the scaling behavior of the 
system and give a value for the exponents T = 1.6 ±0.1 and 
(7 = 0.91 ±0.15. This value for T is surprisingly close to the 
mean-field value predicted for systems of four or more dimen- 
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FIG. 6: The moments (s) and (s 2 ) of the avalanche size distribution 
as a function of (a) drop height h above the apex of the pile for glass, 
stainless steel and zirconium beads and (b) width w across the pile 
for zirconium beads. The lines show an h~ l l 2 (or w -1 ' 2 ) dependence 
for (s) and an h~ 3 ' 2 (or w~ 3 / 2 ) dependence for (s 2 ). The power-law 
behaviors confirm scaling in this system 
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